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■ Abstract. Interstellar absorption in the galactic plane is 
, highly variable from one direction to another. In this pa- 

■ per colour excesses and distances from a new open cluster 
I sample are used to investigate the spatial distribution of 

■ the interstellar extinction. An inverse method (Tarantola 
[ & Valette, 1982) is used to construct the extinction map 
' in the galactic plane below |5| < 10°. The Ay{rJ) dia- 
, grams are compared with those derived from individual 

■ stars (Arenou et al. 1992, Neckel & Klare 1980). An ana- 
. lytic expression for the interstellar extinction as a function 

■ of galactic longitude and distance in the solar neighbor- 
, hood is given. The comparison of the model predictions 

■ with Hipparcos observations in the 4-dimensional space of 
' {V, B — V, Hy, r) shows that our extinction model pro- 
' vides a better fit to the data. However, a new and more 
I detailed extinction model is still lacking. 

Key words: Interstellar extinction ~ solar neighbourhood 

' - Hipparcos - Methods: statistical - open clusters and 

' associations: general 



1. Introduction 

The light extinction caused by the interstellar medium in 
the solar neighborhood alters the apparent brightness of 
the stars and thus their intrinsic properties (absolute mag- 
nitude, intrinsic colour, etc.). For this reason, it is crucial 
to know the effects of the interstellar medium on the ra- 
diation we receive from the stars, and to build models up 
which reasonably predict the amount of extinction in a 
given direction. 

Send offprint requests to: bclien@mizar.am.ub.es 
* Based in part on observations made with the ESA Hippar- 
cos astrometry satellite 



Bahcall & Soneira (1980) have used Sandage's extinc- 
tion law (1972) in their model of the Galaxy, and com- 
pared model predictions with star count observations in 
17 fields. They found that, at | b | > 10°, the Sandage 
absorption law is a good approximation to observations. 
However, at low galactic latitudes, the absorptions are 
patchy and poorly known. In the context of our study on 
stellar kinematics from Hipparcos observations and open 
cluster (Chen et al, 1998), we need a reliable extinction 
model at low galactic latitudes. 

Neckel & Klare (1980) derived extinctions and dis- 
tances for more than 11000 stars and investigated the 
spatial distribution of the interstellar extinction at | b | < 
7.6°. Their results are shown in a series of diagrams. Us- 
ing a larger sample (about 17000 stars) with MK spectral 
types and photoelectric photometry, Arenou et al. (1992) 
(hereafter AGG92) constructed an extinction model in 
which the sky was divided into 199 cells. In each cell, 
a quadratic relation was adopted. This extinction model 
will be the reference in this paper. Hakkila et al. (1997) 
have developed a computerized model of large-scale visual 
interstellar extinction. The code merges several published 
studies (FitzGerald 1968; Neckel & Klare 1980; Berdnikov 
& Pavlovskaya 1991; AGG92; etc.), and can be used for 
making corrections to individual observations and for cor- 
recting statistical sample. With extinctions determined 
from Stromgren photometry and parallaxes from Hippar- 
cos for a sample of 3799 stars, Vergely et al (1997) have 
constructed a three dimension distribution of absorbing 
matter in the solar neighbourhood. However, all these au- 
thors use similar dataset (stars along with their corre- 
sponding colors, magnitudes, and spectral/morphological 
classification types) so that the results are not statistically 
independent (Hakkila et al. 1997). 

In this paper, we are going to use open clusters rather 
than individual stars to investigate the distribution of in- 
terstellar extinction, and compare our results with those 
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found in the literature. The reason for this investigation is 

twofold. First, the advent of CCDs with their sensitivity 
and linearity has allowed to obtain deep photometry of 
open clusters (Carraro & Ortolani, 1994, Kaluzny 1994). 
Second, distances and reddening for open clusters are more 
reliable than those from individual stars, simply because 
one is dealing with a group of stars and therefore the re- 
sults are less sensitive to individual errors. Thus the open 
clusters system provides us with another way to obtain 
informations on A„ as a function of the galactic longitude ^ 

1 and distance r. ■§ 

Pandey & Mahra (1987) firstly used the open clus- B 
ter catalogue of Janes & Adler (1982) to investigate the ^ 
spatial distribution of interstellar extinction within 2 kpc. 
Their results are presented in graphical form. No analyti- 
cal expression for the interstellar extinction has been pre- 
sented. In this paper, we have used a new open cluster 
database and an inverse method (Tarantola & Valette, 
1982) to derive an analytical expression of the interstel- 
lar extinction in the solar vicinity, and compared different 
extinction models with Hipparcos observations. In section 

2 we present the data used in this paper. In section 3 we 
briefly explain the AGG92's model and compare it with 
our open cluster data. In section 4, we construct a new 
extinction model from open cluster system. In section 5 
the predictions from AGG92's extinction model and our 
extinction model are compared with the Hipparcos data. 
The main conclusions are summarized in section 6. 
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2. The Data 

The main source of the data in this investigation is 
the Open Cluster Data Base (Mermilliod 1992). In this 
database, the reddening and distance for each cluster were 
determined from a system of weights corresponding to the 
precision of open cluster data. 412 open clusters with red- 
dening and distance are listed in the database. In addi- 
tion to this data, from the list of open clusters, we have 
checked each cluster which has no information about red- 
dening and distance from ADS database to see whether 
there are new observations published after 1992. Finally, 
we have found 22 such clusters, mostly due to the observa- 
tions from Padova and Boston group (Priel, 1995). Thus, 
the present study is based on the observed distances and 
colour excesses of 434 open clusters. The ratio of total to 
selective absorption, i?=3.1, was adopted to obtain total 
absorption from selective absorption E{B — V). In Figure 
1 we show the distribution of the open clusters in the (l, 
b) plane of galactic longitude and galactic latitude. We 
can see that most open clusters (91%) lie on the galactic 
plane ( | b | < 10° ), and that the distribution of open 
clusters is not homogeneous in galactic longitude (l). In 
the direction between 1= 30° and I = 90°, and between I 
= 145° and I — 220°, fewer open clusters are found. 

In Figure 2 we show the histogram of the distance r. 
In the Open Cluster Data Base (Mermilliod, 1992), the 
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Fig. 1. The distribution of the open clusters in the plane of 
galactic longitude and latitude 



errors of the extinction and distance are not given. We 
suppose that the error aAv can be written : 



f^Aviobs) = AAv{obs)cal + AAv{obs)r 



(1) 



Where AAv{obs) cai is the calibration error approxi- 
mately equal to 0.1 mag and AAv{obs)r the extinction 
error due to the distance error. As a first approximation 
we consider that the distance determination and Av{obs) 
are uncorrected, and that the relative distance error is 
0.25. Then, we have 



<^Av(obs) = 0.1 + 0.25Av{obs) 



(2) 



From Hipparcos parallaxes, Mermilliod et al (1997) 
and van Leeuwen & Hansen Ruiz (1997) claimed that the 
distance of young Pleiades cluster is 10^^ 15% smaller than 
its photometric distance. In order to check whether a sys- 
tematic bias exists on all determinations of open clusters 
from photometric methods, we have compared the open 
cluster distance derived from Hipparcos data (Robichon 
et al. 1997; Mermilliod et al. 1997) with those in Mer- 
milliod's catalogue. Figure 3 shows the distance derived 
from Hipparcos parallaxes (ri) against those in Mermil- 
liod's open cluster database (r2), no systematic trend has 
been found. 
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3. AGG92's extinction model 
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AGG92 have used INCA database (Gomez et al. 1989; 
Turon ct al. 1991) to construct a tridimensional model 
of the galactic interstellar extinction. For each star the 
distance r to the sun and the visual extinction Ay have 
been computed from MK spectral type and photoelectric 
colour index. 
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Fig. 2. Histogram of distance r to the Sun of the open clusters 
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Fig. 4. The observed extinction A„(obs) against that derived 
from AGG92's model A„ (model) for our open cluster sample 



The sky was divided in 199 cells. In each cell, based on 
the work of Neckel & Klare (1980), a quadratic relation 
was adopted: 



a Per 
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Ay = ar + /3r'^ 



(3) 
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Fig. 3. The open cluster distance derived from Hipparcos data 
(rl) against those in Mermilliod's open cluster database (r2) 



For each cell an analytic expression of the interstel- 
lar monochromatic extinction at V-magnitude A„(r,Z,6) 
was obtained. We have used AGG92's model to derive the 
extinction at the position of each open cluster in our sam- 
ple. In Figure 4, we plot the observed extinction j4^(obs) 
against that derived from AGG92's model (model). It 
can be seen that AGG92's model underestimates large ex- 
tinction value. Since AGG92 used only stars brighter than 
V = 13 mag, as they stated, their model cannot be ex- 
trapolated farther than 1 kpc. From Figure 4 we have also 
found that AGG92's results are higher than the observa- 
tions for small distance (r < 500 pc). AGG92 pointed out 
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Extinction structures in the galactic plane. 
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Fig. 5. Extinction structures in the Galactic plane 



that a slight overestimation of Ay for small distance is ex- 
pected in their model, this is consistent with what we have 
noticed from our open cluster sample. In the following sec- 
tions, we use open clusters with | b | < 10° to derive the 
interstellar extinction at low galactic latitudes and con- 
struct a new extinction model in solar neighbourhood. 

4. Using open clusters to derive the interstellar ex- 
tinction at low galactic latitudes 

In order to determine the opacity in each point of the 
space in a circle of 2500 pc around the Sun, we utilized 



an inverse method (Tarantola & Valette, 1982; Valette & 
Vergely, in preparation) which for the first time has been 
applied to detect extinction structures in the Hipparcos 
data (Vergely et al., 1997). 

In this study, we have assumed that the extinction fol- 
lows a linear law : 



Jo 



r)dr 



where Pop{r) is the opacity at the point r in the V-band 
(in mag/pc). Ei is the integrated extinction in the V-band 
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Fig. 6. Extinction structures in the cross section of the galactic plane 
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Fig. 7. Histogram of 5 value {6 = ±^(25£)-i^i^i2^). 



along the line of sight {li, bi, ri) in galactic coordinates. 
The extinction decreases strongly with the galactic lati- 




Av(obs) 

Fig. 8. The extinction derived from our model Av(model) 
against the observed one Av(obs). 



been chosen : 



pop{l, b, r) = Ao{l, b, r) exp ( - 



|r sin(6)| 
ho 



Aq{1, b, r) represents the extinction fluctuations around 
tude. For this reason, an opacity exponential model has an exponential mean extinction law at the point {l,b,r) and 
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Fig. 9. The extinction as a function of distance for each of the 36 cells. The solid lines are the results in this study from the 
inverse method, the dashed lines are the results from Arenou et al. (1992). The crosses (+) are the open cluster sample. 
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Fig. 9. (continued)The extinction as a function of distance for each of the 36 cells. The solid lines are the results in this study 
from the inverse method, the dashed lines are the results from Arenou et al. (1992). The crosses (+) are the open cluster sample. 
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Table 1. The coefficients of the extinction curves in different longitudes 



Range of / 


a I 


Cl2 




04 


(degrees) 


(mag.kpc"'^) 


(mag.kpc"^) 


(mag.kpc""*) 


(mag.kpc"'') 


0-10 


0.752 


0.610 


0.090 


-0.387 


10-20 


0.500 


3.012 


-3.659 


1.378 


20-30 


0.245 


5.468 


-7.721 


3.432 


30-40 


0.156 


5.996 


-9.214 


4.478 


40-50 


0.300 


3.753 


-5.793 


2.893 


50-60 


0.490 


9.364 


-1.023 


0.476 


60-70 


0.740 


-1.825 


3.780 


-1.622 


70-80 


0.943 


-3.715 


6.975 


-2.682 


80-90 


0.747 


-2.954 


6.672 


-3.190 


90-100 


0.598 


-2.604 


7.809 


-4.593 


100-110 


0.442 


-1.745 


7.066 


-4.558 


110-120 


0.235 


0.201 


2.681 


-2.111 


120-130 


0.065 


1.821 


-0.625 


-0.308 


130-140 


-0.103 


3.096 


-2.055 


0.248 


140-150 


-0.248 


4.273 


-3.300 


0.440 


150-160 


-0.189 


4.592 


-4.951 


1.556 


160-170 


-0.033 


4.066 


-5.594 


2.389 


170-180 


0.091 


3.101 


-4.644 


2.158 


180-190 


0.163 


2.070 


-3.064 


1.476 


190-200 


0.217 


1.084 


-1.410 


0.698 


200-210 


0.286 


0.189 


-0.117 


0.392 


210-220 


0.307 


-0.227 


0.434 


0.288 


220-230 


0.259 


-0.160 


0.649 


-0.266 


230-240 


0.201 


0.192 


0.123 


-0.128 


240-250 


0.198 


0.258 


0.160 


-0.136 


250-260 


0.274 


-0.082 


0.729 


-0.207 


260-270 


0.313 


-0.016 


0.187 


0.374 


270-280 


0.288 


0.386 


-0.462 


0.361 


280-290 


0.313 


0.535 


-0.650 


0.395 


290-300 


0.381 


0.398 


0.079 


-0.117 


300-310 


0.280 


1.661 


-1.839 


0.643 


310-320 


0.120 


3.513 


-5.033 


2.260 


320-330 


-0.046 


5.282 


-7.844 


3.611 


330-340 


-0.031 


5.477 


-7.961 


3.546 


340-350 


0.465 


2.076 


-2.138 


0.871 


350-360 


0.961 


-1.226 


3.345 


-1.648 



ho, the characteristic height of the extinction structures. 
The value of /iq is a parameter which is determined dur- 
ing the inverse process. The computed value of h{) is 120 
pc that is consistent with previous studies (e.g. Sharov, 
1964). Since the sample of clusters is limited to a finite 
number of linos of sight that arc relatively spread-out, wo 
have assumed some regular properties of pop{l,b,r) and 
introduced a smoothing length which corresponds to the 
mean intcr-clustcrs distance. The details of the absorption 
clouds, with a length lower than the correlation length, are 
not detected. A 300 pc correlation length gives us global 
fluctuation extinction tendencies. In Fig. 5 and Fig. 6, we 
show the extinction structure in the galactic plane and in 
the cross section of the galactic plane, respectively. These 



figures show that the absorption in the galactic plane is 
highly variable from almost zero to 1.5 mag/kpc, no in- 
formation are available at high galactic latitude. 

Fig. 7 shows the histogram of 5 value {6 = 

— — -), and m ligure 8 we plot the observed 

extinction Av(obs) against that derived from our model 
Av(model). Prom Figure 8, we did not see an overesti- 
mation A^, (model) for small distance (r < 500 pc), which 
seems to be the most significant improvement in our model 
over AGG92's model (see Figure 4). The error distribution 
in Figure 7 is asymmetrical and the negative tail shows 
that some extinction values are overestimated. This could 
be explained by the fact that the model does not take into 
account the fluctuations at the smallest scales and does 
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not represent windows without extinction. Prom Figure 7 

and Figure 8, we can also see that model underestimates 
some extinction values at large Av(obs). Recently, Hakkila 
et al. (1997) found that large-scale extinction studies tend 
to underestimate extinction at distances r > I kpc. This 
is in agreement with our result. For this reason, in the 
following sections, we limit our studies for r < 1000 pc. 

In order to establish an extinction model in the solar 
neighborhood, we have divided the galactic plane into 36 
cells with AZ = 10°. At the center of each cell, the extinc- 
tion curve has been computed. Figure 9 shows the A-^ as 
a function of r for each of the 36 cells. We have compared 
Figure 9 with that of Neckel & Klare (1980) and Pandey 
& Mahra (1987), and found that, although A„(r) changes 
from one cell to another significantly, the scatter in the 
A„ (r) is in agreement with that in the diagrams of Neckel 
& Klare (1980) and Pandey & Mahra (1987). By compar- 
ing our results with those of AGG92, we can see that the 
general trends are similar: strong extinction values in the 
longitude interval 20° — 40° and week extinction values in 
the interval 210° — 240°. However, we have observed the 
presence of a strong extinction in the 350° — 360° interval 
at about 1.5 kpc which is not evident in AGG92's model. 

For each cell, an analytic polynomial expression is used 
to fit the results from the inverse method within r < 1000 
pc : 

E{1) = air + a2r^ + aar^ + a^r'^ (4) 

where E{1) is the extrapolated extinction at the longitude 
I which belongs to the interval [lAl, {i + l)Al] and r is 
the distance in kpc. The coefficients 01,02,03 and 04 are 
given in Table 1. 

In order to provide a formal precision on our extinction 
model, we have made some comparisons with the observa- 
tions for 155 stars from Neckel & Klare (1980). The stars 
selected have a spectral type B, A or F. 0-type stars are re- 
jected because they could have circumstellar extinction. In 
Figure 10, we plot the observed extinction A„(obs) against 
that derived from our model predictions A„ (model). Con- 
sidering the error on extinction the agreement is satisfac- 
tory. We have grouped the Neckel & Klare sample into four 
groups according to A^ (model). In Figure 11, we show the 
histograms of Av(obs) in each group. Since the sample is 
small in each group, the errors {(7^^) have been derived by 
a robust estimator (Morrison & Welsh 1989; Chen 1997a). 
Table 2 shows the results. 

Where (tav is the total extinction error on A^,. As a 
first approximation we neglect the contributions of the 
photometric errors in the Neckle & Klare data, we found 
a Avi'model) = 0.17 + 0.38 Av(model), which gives an ap- 
proximated estimation of the external error in our extinc- 
tion model. 

For |b| > 10° we adopt the Sandage absorption law, 
which is widely used through the literature. Finally, our 
extinction law in solar neighborhood is constructed as fol- 
lows. 



Table 2. The extinction error (ctav) in each group Av(model) 



Av{model) Av{model) ctav number of stars 



0-0.25 


0.15 


0.20 


60 


0.25-0.5 


0.35 


0.27 


31 


0.5-1.0 


0.75 


0.56 


42 


1.0-1.5 


1.21 


0.57 


22 



Ay{r,l,b)=0 \b\>50 

A^ir,l,b) = 0.165(1.192- | tanb |) | cscb \ [1 - exp{-r | 

sinb\h-^)] 10<|6|<50 

A,{r,l,b) =Eil) \b\<10 

Zero-reddening value in Sandage's model for | b | > 
50 has been criticized by de Vaucouleurs & Buta (1983), 
and Lu et al. (1992), who found a polar extinction of 0.1 
to 0.2 mag in A^. Prom the color excesses derived from 
Stromgrc^ii photometry for a complete sample of A3-G0 
stars brighter than B = 11.5 and | b | > 70, Knude (1996) 
found that the North Galactic Pole exhibits a complex 
reddening distribution. However, we found that Sandage's 
model with |b| > 50 is in agreement with the Burstein 
& Heiles maps (1982). Moreover, Hakkila et al. (1997) 
compared several previous results in the directions of the 
Galactic poles (see their Table 3), and found an average 
extinction A„ = 0.1 ± 0.2 mag in both galactic poles. We 
believe that Sandage's model in high galactic latitudes is 
a good approximation to observations. 

5. Comparison with Hipparcos observations 

5.1. Hipparcos observations 

Hipparcos provides an immense quantity of accurate astro- 
metric and photometric data from which many branches 
of astrophysics are benefiting. The primary result is an 
astrometric catalogue of 118218 entries nearly evenly dis- 
tributed over the sky with an astrometric precision in po- 
sition, proper motion and parallax of 1 mas or mas/yr, or 
better for the brightest stars. More details can be found 
in the comprehensive introduction of the catalogue (Fer- 
ryman, 1997, ESA, 1997). Hipparcos observations include 
a survey, which is a complete sample selected according 
to a combined magnitude-spectral class -galactic latitude 
criterion. The selection criterion is V < 7.9 4- 1.1 | sin(b) 
I for stars earlier than or equal to G5, and V < 7.3 -|- 1.1 
I sin(b) I for later spectral type stars. 

In this paper, we have used a complete sample with V 

< 7.3 mag. Early-type stars are known to be associated 
strongly with the Gould belt and the moving groups (Chen 
et al. 1997). The non- uniform distribution of early type 
stars is not included in our Galaxy model. Therefore we 
have excluded the O, B and earlier A- type stars with B-V 

< 0.0. With this selection criterion we are left with 18072 
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Av(rnodel) 



Fig. 10. The extinction derived from our model against the 
observations for 155 stars from Neckel & Klare (1980) 



V > from the Hipparcos 



stars with V < 7.3 mag and B 
catalogue. 

5.2. Galactic model 



A Galaxy structure and kinematic model has been con- 
structed (Chen 1997b). The model includes a thin disk, a 
thick disk, and a halo, which can predict the magnitudes, 
positions, colours, proper motions, radial velocities and 
metallicities according to the selection criteria used in the 
observations. Since disk stars are dominating in Hipparcos 
survey, we simply describe the main characteristics of the 
disk population in our Galaxy model. Any further detail 
of the model can be found in Chen (1997b). 

The stellar density laws used for the disk population 
are exponential. For the distribution of stars perpendic- 
ular to the plane of the Galaxy, which is known to vary 
with luminosity, we adopted the relations derived by Bah- 



call & Soneira (1980). The luminosity function for the disk 

stars given by Wiclcn ct al. (1983) has been used in our 
model. The variation of the velocity dispersion in the so- 
lar neighborhood with spectral type and luminosity class 
has been adopted from Dclhaye (1965) and Ratnatunga et 
al. (1989). We have adopted a velocity dispersion gradi- 
ents from the Galactic plane suggested by Fuchs & Wielen 
(1987). Mendez & van Altena (1996) have derived the dis- 
tance of the Sun from the Galactic plane, Z = 2 ± 34 pc 
and Z = -8 ± 19 pc, from two fields. But, recently they 
(Mendez & van Altena, 1998) found a distance of the Sun 
from the symmetry plane of the Galaxy of Z = 27 ± 3 pc. 
Using a large sample of OB stars within 10 degrees of the 
galactic plane. Reed (1997) has investigated the Sun's dis- 
placement from the galactic plane and found Z = 10 12 
pc. We found that this parameter is not well determined, 
for this moment, we keep Z = in our model. 

In order to make comparisons with the Hipparcos ob- 
servations, an area integration was performed by adding 
up the contributions of many regions nearly uniformly dis- 
tributed over the sky. A grid with spacings of 5° in lati- 
tude and 15° in longitude, divided the sky into 562 regions 
each of 75 dcg^. Longitude spacings at high latitude are 
increased from the value selected at the Galactic equator 
to keep the areas in each cell approximately constant. 

5.3. Comparisons of the model with a selected sample 
from Hipparcos data 

Domingo (1998) has compiled a sample of B, A and F-type 
stars from Hipparcos observations with uvbyH^^ photom- 
etry. Distances have been derived from Hipparcos paral- 
laxes, and the Stromgren photometric data come from the 
Hauck and Mcrmilliod (1990) compilation and the new ob- 
servations performed by Barcelona group (Figueras et al., 
1991, Jordi et al., 1996). In order to derive reliable phys- 
ical parameters, we have not included in our sample any 
stars known to be or suspected of being variable, spectro- 
scopic binary, known as peculiar (Am, Ap, 6 Del, ...) and 
those with variable radial velocity. Stars belonging to dou- 
ble or multiple systems for which only joint photometry is 
available have also been rejected. 

The absorption has been derived from Stromgren pho- 
tometry (Domingo, 1998). We have selected stars in our 
sample with A„ > and r < 1000 pc. Our sample includes 
1492 stars and has an average distance f = 177 pc. In 
Figure 12, we have compared AGG92 and our extinction 
model with the observations. The average extinction de- 
rived from Stromgren photometry in the sample is 0.121 
± 0.003 mag, and the average extinctions derived from 
AGG92 and our extinction model are 0.203 ± 0.005 mag 
and 0.127 ± 0.003 mag, respectively. It can be seen that 
AGG92 model provides a systematically larger extinction, 
and our model gives a better fit to observations. This is in 
good agreement with our result from open cluster dataset. 
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Fig. 11. The histogram of Av(obs) 



5.4- Comparisons of the model with Hipparcos data from 
whole sky 

The predictions from our and AGG92 extinction models 

have been compared with the Hipparcos observations in 
Figure 13. The reduced proper motion Hy is defined ac- 
cording to Luyten (1922): 

Hy=V + 5\ogfi + 5 (5) 

where /i is the total proper motion expressed in arcsec 
yr~^. With our extinction law, the model predicts 18486 
stars, 2.3% more than that from the Hipparcos observa- 
tions, while with the Arcnou ct al. (1992) extinction law, 
the model predicts 17229 stars, 4.7% less than that from 
real observations. The main features of the observed dis- 
tribution, for example, the two peaks near B-V = and 



1.0 mag, are present in the models and have the correct 
amplitude. In order to compare different extinction laws 
with observations in the 4-dimensional space of {V, B — V, 
Hy, r), we use a cluster analysis algorithm (Chen 1996) 
to choose the best fit model. The Hipparcos data, made 
of magnitude {V), colour (B — V), reduced proper motion 
(Hy) and distance (r), is merged with the simulated sam- 
ple. We carry out cluster analysis for this merged sample, 
then, we separate the observed data from the simulated 
one and compare the distribution between the model pre- 
dicted stars and the real observed stars. The content of 
clusters can be interpreted from the model predictions. 
Table 3 gives the main physical parameters of each clus- 
ter from OTU- model predictions, and Figure 14 shows the 
positions of each cluster in H-R diagram. 
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Fig. 12. Comparison between the observed A„ distributions and the model predictions, a): AGG92 model has been adopted, 
b): Our model has been adopted 



Table 3. The means and standard deviations (between brackets) of observables and physical parameters in each cluster together 
with the cluster size. The physical parameters (My, au:cv-'crw) are derived by the corresponding simulated stars 





V 


M„ 


B-V 


r 


cru'crv'-O'w 


N 




mag 


mag 


mag 


pc 


km/s 


stars 


Cluster 1 


7.0 


1.3 


0.3 


154 


21:14:13 


3948 




(0.21) 


(1.64) 


(0.23) 


(96.2) 






Cluster 2 


4.0 





1.0 


80.4 


34:21:17 


782 




(0.89) 


(1.83) 


(0.54) 


(66.1) 






Cluster 3 


6.4 


2.3 


0.5 


84.7 


27:19:16 


3049 




(0.33) 


(1.73) 


(0.29) 


(58.3) 






Cluster 4 


6.9 


0.67 


1.15 


182 


37:28:20 


3239 




(0.29) 


(1.27) 


(0.34) 


(100.8) 






Cluster 5 


5.3 


1.7 


0.4 


62 


25:17:14 


1535 




(0.47) 


(1.67) 


(0.35) 


(44.7) 






Cluster 6 


5.8 


-0.4 


1.5 


178 


37:24:19 


2357 




(0.44) 


(1.20) 


(0.29) 


(84.0) 






Cluster 7 


6.9 


-1.4 


1.6 


420 


42:29:23 


3162 




(0.35) 


(0.98) 


(0.39) 


(163.4) 
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Fig. 13. Predicted vs. Observed {circles) distributions for whole sky. The solid line indicates our model predictions, while the 
dotted line indicates the prediction from AGG92's model. 



The statistics is used to test the capability of mod- 
els to represent the data and choose the best fit to the 
observation: 

i=l model 

where iV*j,^ and N^^^^i arc the number of stars in clus- 
ter i in observed data and simulated data, respectively, k 
is the total number of clusters, we found that fc = 7 is 
good enough for this comparison (Chen 1996). 



Obviously, if the model is a good representation of the 

observations, then it should produce frequencies compa- 
rable with the observed ones in each cluster. A large value 
of indicates that the null hypothesis is rather unlikely. 
From the sample, we have derived = 28 and 43 for our 
extinction model and AGG92's model. The method (called 
as chi-by eye) described above has been tested by Monte 
Carlo simulations and used to select the best-fitting model 
(Chen 1996). This method has also been used by Ojha et 
al(1996) to constrain galactic structure parameters. Our 
results show that our extinction model gives a better fit to 
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Hipparcos observations. We should point out that in our 

and AGG92 extinction model, the probabilities that the 
model reflects the reality are ruled out at 3.4 and 4.8 sigma 
Gaussian equivalent level. This is due to the fact that the 
statistics of the errors in the data is not a Poisson statis- 
tics, because of some systematic errors in the photometry 
and astrometry in our Galaxy model, such as disk colour- 
magnitude diagram, absolute magnitude calibration for 
giants, velocity distribution and density law of the galac- 
tic disk stars, and the little fluctuation of the interstellar 
clouds as well. In this paper, for a given Galaxy model 
(Chen 1997b), we compare model predictions from differ- 
ent extinction laws with observations, which can give us 
some information about the uncertainty in the simulated 
sample due to the extinction. In a forthcoming paper, we 
plan to improve galactic structure and kinematical param- 
eters from Tycho catalogue with the help of our Galaxy 
model and extinction law. 



80 



B 40 - 




-80 



Kfyr 




Isochrones for Z=0. 02 
-0.5 0.0 0.5 



2.0 



B-V 



Fig. 14. The means and standard deviations of 7 clusters. 
Isochrones from Schaller et al. (1992) for solar metal content 
have been overplotted. Considering that our sample has an 
average extinction A,, = 0.2 to 0.3 mag, B-V for Isochrones 
has been shifted 0.1 mag from (B-V)o. 



5.5. Comparisons of the model with Hipparcos data at low 
galactic latitudes 

In this section, we compare the model predictions with 
Hipparcos observations at low galactic latitudes (| b | < 
10°). Since the stellar density is about 5.5 stars per square 
degree in the galactic plane from Hipparcos observations. 



10 20 30 
area number 



40 



Fig. 15. The difference 5 = N{obs) — N(model) as a function 
of the cell number (from 1 to 36). The solid line indicates our 
of the cell number (from 1 to 36). The solid line indicates our 
model predictions, while the dotted line indicates the predic- 
tion from AGG92's model 



comparisons must involve relatively large areas in the sky. 
We have used 36 areas, each of them covers from -10° to 
10° in galactic latitude and 10° in galactic longitude. 

AGG92's extinction law and the law presented in this 
study have been used for comparisons with the observa- 
tions. In Table 4, we show the number of stars in each cell 
derived from Hipparcos observation (A^o) and predicted 
by AGG92's extinction law (A''i) and our extinction law 
(A^2)- We have used the statistics (eq. 6) to compare 
the model predictions and observations. With our extinc- 
tion law, we found = 62.3 and for AGG92's model, x^ 
= 108.7. The model predictions are at 2.6 and 5.7 sigmas 
from the data. 

In Figure 15, we plot the difference 5 = N{obs) — 
N {model) as a function of the cell number (from 1 to 36). 
From x^ and Figure 15, we can see that our extinction 
law provides a better fit to Hipparcos observations than 
AGG92's one, which predicts less stars than observed. On 
the another hand, we found that both models predict too 
less stars between ^ = 50 to 130. From Table 4, we can 
see there are 152 stars between I = 280 to 290, and only 
80 stars between I = 150 to 160. Neither our extinction 
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Table 4. The number of stars in each cell derived from Hipparcos observation (Nq) and predicted by AGG92's extinction law 
(Ni ) and our extinction law (N2 )■ 



Range of I No Ni N2 Range of / A^o Ni N2 

(degrees) (stars) (stars) (stars) (degrees) (stars) (stars) (stars) 



0-10 


97 


111 


115 


180-190 


103 


107 


115 


10-20 


104 


87 


109 


190-200 


119 


114 


121 


20-30 


104 


87 


105 


200-210 


135 


116 


131 


30-40 


91 


106 


112 


210-220 


107 


108 


114 


40-50 


121 


97 


120 


220-230 


114 


105 


112 


50-60 


131 


103 


115 


230-240 


125 


107 


126 


60-70 


142 


119 


121 


240-250 


95 


101 


127 


70-80 


133 


115 


123 


250-260 


137 


117 


128 


80-90 


136 


116 


120 


260-270 


132 


113 


121 


90-100 


116 


109 


117 


270-280 


112 


114 


119 


100-110 


151 


118 


123 


280-290 


152 


113 


121 


110-120 


135 


112 


126 


290-300 


150 


115 


122 


120-130 


133 


116 


121 


300-310 


132 


112 


120 


130-140 


120 


108 


117 


310-320 


124 


109 


123 


140-150 


127 


113 


116 


320-330 


126 


107 


124 


150-160 


80 


101 


103 


330-340 


101 


105 


118 


160-170 


100 


103 


107 


340-350 


109 


107 


123 


170-180 


123 


113 


119 


350-360 


117 


105 


119 



law nor AGG92's model can model this large change. A 
more detailed absorption law is needed for modeling the 
extinction substructure in the galactic plane. 

6. Discussions and Conclusions 

At low galactic latitudes (| b | < 10°), the absorptions 
are patchy and not well known. In this paper colour ex- 
cesses and distances of open clusters compiled by Mer- 
milliod (1992) have been used to study the interstellar 
extinction. We have compared our extinction model with 
AGG92 and Hipparcos observations. The main results of 
this paper can be summarized as follows. 

1. An inverse method (Tarantola & Valette, 1982) was 
used to construct the extinction map in the galactic plane. 
We found the lowest absorption between Z = 210° to 240°, 
and the highest absorption between I = 20° to 40°. 

2. By using three different samples, we found that 
AGG92's results are higher than the observations for small 
distances. 

3. We have constructed an extinction law by combining 
Sandage extinction law with | b | > 10° and our extinction 
law in the galactic plane. 

4. AGG92 and our extinction models have been com- 
pared with a well selected sample from Hipparcos ob- 
servations with Stromgren photometry. We found that 
AGG92's model provides a slight overestimation of A„ . 

5. Galaxy model predictions with our and AGG92's 
extinction laws have been compared with the Hipparcos 
observations for the whole sky and at low galactic lat- 
itudes, respectively. The main characteristics in the ob- 



served distributions (for example, two peaks in B — V) 
are well predicted by both models. The statistics, used 
to choose; the best fit to Hipparcos observations, shows 
that our extinction model provides a better fit to data. 

6. Interstellar absorption in the galactic plane is highly 
variable from one direction to another. We believe that 
a more detailed absorption law, which can represent the 
extinction at smaller scales, is needed for modeling the 
extinction substructure in the galactic plane. 

These results are very useful for us to investigate the 
stellar kinematics in the solar neighborhood from Hippar- 
cos and Tycho data in the near future. But a new and more 
detailed extinction model from Hipparcos observations is 
still missing. 
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